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Abstract 

This paper presents several matrix elements for processes that may con- 
tribute to e + e~ — ► 7r + 7r - 7r°7T derived using the isobar model and the Lorentz 
invariant amplitude method. The formulas may be used to measure the elec- 
tromagnetic form factors of the p meson in experimental data. Quantum 
Chromodynamics (QCD) calculations predict the electromagnetic form fac- 
tors of the p meson which can be compared to the measured values. The 
formulas in this paper may be used to study the a\ meson if e + e~ — ► a\n 

contributes to e + e~ — ► 7r + 7r~7r 7r°. 
23.23.+X, 56.65.Dy 
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I. INTRODUCTION 



The isobar model (described below) is a phenomenological theory of intermediate energy 
(a few GeV) strong interactions ]!]-§]. Phenomenology is "a descriptive or classificatory 
account of the phenomena of a given body of knowledge, without any further attempt at 
explanation" ||. In this sense, the isobar model is phenomenology. Quantum Chromody- 
namics (QCD) hopefully provides a fundamental explanation of the interactions described 
by the isobar model. The isobar model cannot be derived at present from QCD. It presum- 
ably is not a fundamental explanation of strong interactions, but provides a description of 
many phenomena in strong interactions. 
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This paper derives formulas for the isobar model processes that may contribute to the 
reaction e + e~ — ► 7r + 7r _ 7r 7r . These formulas can be compared to e + e~ — ► 7r + 7r _ 7r 7r data 
from experiments such as CLEO, CLEO II, the Beijing Electron Spectrometer (BES), Mark 
III , and DM2. 

The formulas can be used to measure the electromagnetic form factors of the p meson in 
experimental data on e + e~ — > 7r + 7r~7r°7r°. There are at least two predictions of the electro- 
magnetic form factors of the p based on QCD arguments to compare to the experimental 



measurements flO Hl^l . The isobar model formulas can be used to determine the properties 
of the a\ resonance if e + e~ — > ayn is a significant contributing process to e + e _ — ► 7r + 7r~7r 7r°. 
A more detailed discussion of these results can be found in the thesis on which this paper is 
based Hl4 . 



II. THE ISOBAR MODEL 

The isobar model assumes that particle production and decay proceeds via resonances. 
For example, in the process 7p — > ir°p, pion photoproduction, the isobar model assumes the 
matrix element is of the general form 

M = £ AiMEfa, e p )— 1 (1) 

where m ; is the mass of resonance % and r, is the width of the resonance. This is the 
Breit-Wigner formula which originated in studies of atomic and nuclear physics. The isobar 
model uses several variants of the Breit-Wigner formula. The Breit-Wigner formulas are the 
non-relativistic Breit-Wigner formula with constant mass m and width T , the relativistic 
Breit-Wigner formula with constant mass and width, and Breit-Wigner formulas where the 
width r(m) is a function of the mass. 

If the resonance is narrow, there is little difference between a Breit-Wigner formula with 
a constant width Tq = T(mo) and a Breit-Wigner formula with a width that is a function of 
the mass. If the resonance is broad, the Breit-Wigner formula with mass-dependent width 
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T(m) may be significantly different from a Breit-Wigner with constant width. Additional 
assumptions are required to specify the mass- dependent width T(m). 

MEi{e 1 , e p ) is a Lorentz invariant matrix element incorporating the spin and momentum 
dependence of the interaction. Ai in Equation [I] could be a complex function of the mass 
m2 = G°7 + Pp) 2 but it is generally a complex constant. This is certainly a reasonable 
approximation if the resonance i is narrow. If the resonance is broad, this approximation 
may be invalid. 

Pion photoproduction is presumed to proceed by 

IP ~* Resonance I — > 7i°p 



In the isobar model, the S matrix [|15|,[nj is the sum of products of Breit-Wigner fac- 
tors, each corresponding to a resonance, and complex matrix elements incorporating spin- 
dependence. Physically, this means the interaction proceeds by sequential two-body inter- 
actions. For example, A + B + C —>■ D + E + F might proceed by A + B ^ D + X followed 
by C + X — > E + F. X is a resonance and appears as a term 

1 

(Pa +Pb- Pc) 2 ~M 2 X + iY x M x 
The process might proceed by A + C ^ D + Y and B + Y^E + Fas well, where Y is a 
different resonance. 

For particle decay, the isobar model is equivalent to sequential two-body decay of res- 
onances. For example, the decay A — > B + C + D is actually A — > X + D followed by 
X — > B + C. The intermediary resonance X appears as a Breit-Wigner term with mass nix 
and width Tx in the S matrix. 

The Breit-Wigner formula is an assumption of the isobar model. This formula originated 
of atomic and nuclear physics. Narrow peaks in cross-sections occur in atomic and nuclear 
physics. These are typically interpreted as bound states in a central potential. The Breit- 
Wigner formula was derived in this context by applying the mathematics of a forced damped 
oscillator such as a tuning fork to atomic and nuclear systems. 



When resonances were discovered in pion-nucleon scattering, an analogy was made to 
the much narrower states observed in atomic and nuclear physics. As part of this analogy 
the Breit-Wigner formula was applied to low-energy hadronic interactions. 

The isobar model incorporates a Lorentz invariant matrix element to represent the angu- 
lar distributions in the interaction due to the spin of the interacting particles. The evaluation 
of these matrix elements comprises the bulk of this paper. 

The isobar model assumes that the resonances have the same properties — the same 
mass, the same width, the same charge, the same quantum numbers — in all processes. 
The mass and width refer to the parameters of the Breit-Wigner formula for the resonance. 
The isobar model assumes a fundamental spectrum of resonances with unique properties 
that describe all low-energy hadronic interactions. The same resonances occur in e + e~ 
annihilation, pion-nucleon scattering, and other processes. In practice, this spectrum is 
finite because only enough energy is present to create the lower energy resonances. 

Generally, isobar models assume that only resonances below a certain cut-off angular 
momentum contribute to low-energy nuclear processes. Resonances with a spin above 2 are 
usually ignored in isobar models. The higher spin resonances presumably have larger masses. 
This makes intuitive sense if the resonances are excited states of quarks in a potential. The 
energy to create the high spin resonances is not available. In practice, radial excitations are 
frequently ignored for the same reason. 

In summary, the isobar model assumes: 

• The S matrix should be a Lorentz invariant function of initial and final state 4-momenta 
and spins. 

• The S matrix can be represented as the sum of products of Breit-Wigners which 
represent known resonances. Physically, this can be thought of as representing the 
process as a sequence of two-body interactions. 

• Resonances are represented by the Breit-Wigner formula. Frequently, the relativistic 
Breit-Wigner formula with constant mass and width is used. 
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• The resonances are the same — have the same mass, width, charge, and other quantum 
numbers — in all processes: e + e _ annihilation, pion-nucleon scattering, kaon-nucleon 
scattering, pion photoproduction, pion electroproduction, and so forth. 

• Each two-body process has an associated Lorentz invariant matrix element that rep- 
resents the angular distributions in the two-body process due to the spins of the inter- 
acting particles. The rules for constructing these matrix elements are described below. 
The helicity formalism is an alternate scheme for constructing these matrix elements 
that is not used in this paper. 

• Each two-body process has an associated amplitude which is treated as a complex 
constant. For spinless particles, the most general Lorentz invariant form of this ampli- 
tude for a two-body process would be an arbitrary function of the two-body invariant 
mass in the two-body process. A complex constant is a reasonable approximation for 
narrow resonances. For broad resonances, it constitutes an assumption. 

In practice, only resonances with less than a cut-off spin J are used and radial excitations 
are frequently ignored. 

III. LORENTZ INVARIANT AMPLITUDE METHOD 

Lorentz invariant matrix elements are constructed by describing a decay as sequential 
two-body or three-body decay vertices. In e + e~ — > n + n~n°n° , p — > nn is a two-body vertex 
and uj — > nnn is a three-body vertex. The total matrix element is the product of the matrix 
elements for each vertex. These matrix elements are constructed from the 4-momenta and 
polarization 4-vectors of the in-going and out-going states. 

In this paper, the following notation will be used. A 4-vector is represented as follows: 

a^ = a= (a x ,a y ,a z ,a t ) (2) 

and 
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a 11 = (-a x , -a y , -a z , a t ) 



(3) 



A 4-momentum is: 



Pv=P= (Px,Py,Pz,E) 



(4) 



where E is the energy. 

This analysis uses two kinds of 4- vectors. 4-momenta change sign under the parity trans- 
formation. Polarization 4-vectors, which are usually represented by e, do not change sign 
under the parity transformation. The amplitudes formed by the Lorentz invariant amplitude 
formalism must have definite parity properties, determined by the intrinsic parities of the 
particles as described below. 

The Minkowski metric g^ v is defined as: 

/-l 0\ 
0-100 
0-10 

V o o 01/ 

Then, the dot product of two 4-vectors can be represented as: 



9yv 



(5) 



(6) 



= 5^b u 



= (a ■ b) 



(a-b) = a t b t - a x b x - a y b y - a z b z 



(7) 
(8) 
(9) 



#f is the Kronecker delta. 



Notice, in this notation: 



p-p = M 2 
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(10) 



With the other possible choice of metric, 

p-p=-M 2 (11) 

The polarization 4-vectors e(A) take the following form in this notation: 

For a massive spin one (1) particle with helicity A with momentum p along the x axis, 



e(A = ±1) = =p(0, 1, ±i, 0)/V2 (12) 

e(A = 0) = (p,0,0,E)/M (13) 
where M is mass of particle, or 

e(X = 0) = (p,0,0,E)/\p\ (14) 
For narrow resonances, Equation ^ and Equation [TJ] are equivalent. However for Equa- 



tion [14], note that if the massive vector is broad, then p ■ p ^ M in general. Equation [14] is 



the correct formula to use in general. Equation |13| is an approximation. These polarization 
4-vectors are orthonormal. They satisfy: 

e(A i )e(A i )* = 8 i3 (15) 

Notice that one dots with complex conjugate to get this. 
For a massless vector particle (such as a photon): 

e-p = e/ = (16) 

The Lorentz invariant amplitudes are also built using the Levi-Civita tensor e^ pcr . The 
Levi-Civita tensor is defined by 

-1, if a(5^8 is an even permutation of xyzt; 
tap-yS = \ +1, if 0(3^5 is an odd permutation of xyzt; 
0, otherwise 



and 



,a/37<5 



The Levi-Civita tensor changes sign under the parity transformation. It has some other 
properties which will be discussed in the section below on the ujtt model. 
The Levi-Civita tensor enters into the amplitudes in forms 



The Levi-Civita tensor is a determinant or a signed volume of a box defined by four 
4-vectors. If a 4-vector is repeated in the above form (e.g. A = B) then the form is 0. 

The matrix elements are built out of the 4-momenta, polarization 4-vectors, Breit- 
Wigners, and the Levi-Civita tensor e^ upa . The polarization of the intermediate states must 
be summed over using the completeness relation. 

PflPu 



(17) 



where e is the polarization 4-vector of the intermediate state, p^ is the 4-momentum 
of the intermediate state, and M is the pole-mass of the intermediate state, is the 
Minkowski tensor. For a narrow resonance, p 2 can be replaced by M 2 in the completeness 
relation. This analysis involves the p which is broad. 

In practice, the completeness relation appears in forms of the type 

(A-P S )(P S -B)- 



J2(A ■ e(X))(B ■ e(X)) = -(A ■ B) 



P 2 



(18) 



where A and B are 4-vectors, either four momenta or polarization 4-vectors. e(A) is the 
polarization of an intermediate state which is summed over. P s is the momentum of the 
intermediate state s. 

The completeness relation also appears in forms of the type 



J2 (A ■ e a (X a ))(e a (X a ) ■ e b (X b ))(e b (X b ) ■ B) 



i Pb-B 
i — ht- 



-(A ■ B) + {A Pb) J"- B) 



(19) 



The Lorentz invariant amplitudes usually reduce to formula constructed of 4-vector prod- 



ucts using the expansions in Equations [Tqand IS above. This allows one to drop the Greek 



indices and provides a form for the amplitudes which is easy to compute with a 4- vector dot 
product function when implementing the model on a computer. 

The specific rules for forming amplitudes (or matrix elements) are: 

• The amplitudes must be formed by contracting combinations of the available 4- 
momenta and polarization 4- vectors. The amplitudes must be Lorentz invariant. 4- 
vectors are not Lorentz invariant, but the contractions are Lorentz invariants. All of 
the polarization 4- vectors must be used to form the amplitudes. The 4-momenta may 
or may not be used. 

• Different amplitudes for the same process must be linearly independent to within 
factors which are rotationally invariant. These factors are envisioned as slowly varying 
functions of invariant masses. 

• The amplitudes must be either scalars or pseudoscalars depending on the required 
parity properties. The product of the intrinsic parities of the decaying particle and 
the decay products determines the required parity transformation properties of the 
amplitude. 

• Photons must satisfy e^p^ = 0. e is the polarization 4-vector of the photon and p is 
the 4-momentum of the photon. 

• Intermediate states must be summed over helicities using the completeness relation 
above. 

• Intermediate states should include a mass dependent phase. This is the Breit-Wigner 
form. 

• If the process being modeled is an eigenstate of C(G) then the amplitude must be 
appropriately symmetrized so that the final state it describes is also an eigenstate of 
C{G). C is charge conjugation which transforms a particle into its anti-particle. G 
refers to G parity. 
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A. Mass Dependent Widths 



The Particle Data Book and other sources only give a single width at the pole mass 
T = r(m ). As mentioned above, the Breit-Wigner form contains a mass dependent width. 
For a particle d which decays to iV particles, the width is determined by 

r(tZ -+ 1, . . . , n) = ^ J dLIPS(m 2 , P Xl . . . , P n )\ME\ 2 (20) 

where ME is the matrix element of the decay. dLIPS refers to integration of Lorentz 
Invariant Phase Space. 

This paper is concerned with two decays: p — > tttt and uo — > tttttt. The volume of phase 
space changes as the invariant mass of the decaying particle d changes. Thus, for this reason 
alone, the width is mass dependent. The isobar model needs formulas that relate T(m) to 
the width r and the mass m of the intermediate state. 

For massive spin one (vector) particles decaying to two particles, the formula 



m p cm 3 1 + (Rp. 



r ( m ) = r o— ) T—TRZ v ( 2l ) 

can be used. It is used for the p resonances. However, it is inappropriate for the u which 
decays to three pions. In this instance, the T(m) can be calculated by explicitly integrating 
the matrix element squared using the formula for the mass dependent width above. 



i2 



In this, J cm0 y2 are the Blatt-Weisskopf barrier penetration factors |17| . The Blatt- 



l + (Rp C m) 2 

Weisskopf barrier penetration factors represent mathematically that the particles are not 
point particles. R is an estimate of the radius of a strongly interacting particle. Notice that 
if R = then the penetration factors are identically 1. The width T is the integral of the 



square of the matrix element of the decay. The Blatt-Weisskopf factor J \+^ cm0 )2 appears 
in the matrix element. This relationship between the width and the matrix element is why 
the square of the Blatt-Weisskopf barrier penetration factor appears in the formula for the 
width. p cm is the momentum of the two-body decay in the resonance rest frame. The value 
of R is not definite. The value 1 fermi or 5 GeV~ x is usually used. 
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IV. LORENTZ INVARIANT AMPLITUDE MODELS FOR E + E~ -> vr+vr 7r°7r° 

This section contains the main results of this paper, Lorentz Invariant Amplitudes for sev- 
eral processes that may contribute to e + e~ — > 7r + 7r~7r°7r°. This paper envisions the process 
e + e~ — > 7r + 7r~7r°7r° as the superposition of several modes listed in Table |I[ All of the modes 
except for the last are represented mathematically as Lorentz invariant complex amplitudes 
that interfere. 

The last mode in Table | is non-resonant four pions, usually production of 7T + 7r~7r 7r° 
according to a four-body phase space distribution added incoherently to the other modes. 
Technically, the isobar model precludes this last mode. It is incorporated in isobar model 
analyses as a substitute for unmodelled modes, detector simulation errors, and so forth when 
performing fits to experimental data. 

A. The e+e" -> p + p~ Mode 

This is V — > VV process where V indicates a vector particle. There are three polarization 
4- vectors e 7 «, e p +, and e p -. e 7 » is the polarization of the virtual photon formed when the 
electron and positron annihilate. There are two independent 4-momenta P and Q. P = 
(0, 0, 0, E cm ) = P p + + P p - and Q = P p + — P p - . There are two neutral pions 7r°. These are 
refered to in this analysis as 7r° for the more energetic pion and for the less energetic 
pion. There are two ways to form the 4-momenta of the two p mesons: P p + = P n + + P n o 
and P p - = P n - + P n o or P p + = P n + + P n o and P p - = P^- + P n o. It will be necessary to add 
the matrix elements for these two cases to insure that the final amplitude has the correct 
Bose symmetry. As a first step, the decay of the p's is ignored. 

There are three independent amplitudes that can be constructed according to the rules 
of the Lorentz invariant formalism. These are: 

ME 1 = (e r -Q)(e p+ -e p -) (22) 



ME 2 = (e r -Q)(e p+ -Q)(e p --Q) 
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(23) 



ME 3 = (e r • e p+ )(e p - • Q) + (e r ■ t p -){t p+ ■ Q) (24) 

Note that these are the matrix elements at the 7* — > p + p~ vertex only. The full formulas 
[see below) must include the p — ► mi decay. All of these formulas are parity — and charge 



conjugation — eigenstates. The ( p + p~) final state is a charge conjugation eigenstate fll8 
Kramer and Walsh's paper quotes these formulas using a different notation. Their paper 
does not consider the p decay. 

Thus, the total matrix element would be 



ME{e + e- -> p + p~) = A\ME\ + A 2 ME 2 + A 3 ME 3 (25) 

where Ai are complex parameters. The parameters correspond to the electric, magnetic, 
and electric quadrupole form factors of the p. The relations between the form factors Fe, 
Fm, and Fq and the A^ amplitudes is 

A 1 = F E + ^F Q (26) 
A 2 = (2(1 - r / ))- 1 [(l - 2 -r,)F Q - F E + F M ] (27) 
A 3 = F M (28) 



where 



p being (0, 0, 0, E cm ) 



B. Incorporating p Decay in Amplitudes 

So far, the propagation and decay of the p's has been neglected. This adds a Breit- 
Wigner factor for each p and a factor e p ■ Q p . Q p ± = P w ± — P^q where the pions are the 
decay products of the p meson. In addition, the polarization of the intermediate p must be 
summed over using the completeness relation. 
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Thus, the matrix elements become 



ME 1 = EA p+ A p _(e 7 * • Q)(t P + ■ e P -)( e P+ ' Q P +)( e P - ' Q P +) 

1 1 



(30) 



( m p+ ~ m l ~ im P + r P + ) ( m p - ~ m l ~ im p- T P -) 
Now perform the summation over the polarizations of the p + and p~ . Recall the relations 

(A.P)(B-P) 



J2(A-e(X))(B.e(X)) = -(A.B) + 

x 



M 2 



J2(A-e a )(e a -e b )(e b - B) 



-(a ■ b) + fA -y 



^(\-(p a -B) + ^- p i^ }) 



Using these relations, the full matrix elements with the p decay become: 



p 



(Qp+- p P +) , 

m 1 



-(P P+ -Q P -) + (P ^- )iP ^- ) 



pT" 

P 



) 



(m 2 -m 2 -im p+ F p+ ) (m^_ -mg -im p _ T p _ ) 



(31) 



(32) 



(33) 



M£ 2 = (e r • Q)([-(Q • Q P+ ) + (Q ' Pp+) ^/ + ' Qp+) 



(Q-fp-)(fp--Q P -) 



P 



) 



ME* 



i i 



( S ,.P p+ )(P p+ -Q p+ ) 



+ 



— (e 7 » • (5/9-) 





p p+ 




-)(Pp-Qp-)] 




p2_ J 








p2 

P 




f)( p p+-Q P +)V 



( m p+ - m o- im P + r P + ) (™ p _ -™d -* m P - r P - ) 



(34) 



(35) 



using these relations, the final expression for the amplitude for the p + p mode(s) becomes 
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+ R 2 e l ^(ME 2 (n° 1 ^) + ME 2 (n 2 n° 1 )) 

+ Rse^iME^n ,) + ME 3 (n°n°)) (36) 

Mi?j(7r°7r2) and MEiij^^i) refer to the two choices of 7r° to form the 4-vectors in order 
to satisfy Bose statistics under interchange of the two identical neutral pions. Defining 7r° 
as the more energetic neutral pion (e.g. E n o > E w o) in the center of mass of the four pions 
frame of reference provides a well-defined and unique way to label the two neutral pions in 
real data. 

C. Simplification of p + p Amplitudes 

The preceding amplitudes are complicated. However, only a few terms dominate the 
amplitudes. The factors (Q p ± ■ P p ±) in the formulas above are the reason for this. 



(Q p ± ■ P p ±) = (iVt - P T o) ■ (iVt + P ff0 ) (37) 



This becomes 



(Q p ± ■ P p ±) = (P n ± ■ P n ±) - (P w o ■ P n o) (38) 



which becomes 



(Q p ± ■ P p ±) = ml ± - ml (39) 

If the masses of the charged and neutral pions were the same, this term would be exactly 
zero. In this case the p + p~ amplitudes reduce to 

ME X = -{e T -Q){Q P+ -Q P -) 

( m l+ - m f- im P + r P +) ( m2 p - - m o- im p- r P - ) V ; 



ME 2 = (e r • Q) (-((Q ■ Q p+ ) + (Q ■ Q p -)) 

i i 

( m l+ - m i~ im p+ r P +) ( m p_ - m i- im P - r P - ) 



(41) 
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ME 3 = - ((e r • Q p+ )(Q -Q p -) + (e 7 * • Q p -)(Q ■ Q p+ )) 

i i ("42) 

( m l+ -■"■>-o- im p+ r P + ) ( m2 p - ~ m o - im p- r P -) V ; 

Although the masses of the charged and neutral pions are slightly different, the terms 
above dominate the p + p~ amplitudes. Another consequence is that the p + p~ amplitudes 
are largely insensitive to the choice of or in the completeness relation. 

D. The e+e -> aivr Mode 

This process consists of 



where 



and 



e + e — > Qq~7r^~ 



af - p ± vr° 



P ± - 7T ± 7T 



At the e+e — > 7* — > ai7r step the relevant 4- vectors are e 7 , e Ql , P = P Ql + Pr and 
Q = P ai — P,-. There are two independent amplitudes: 

MP/ +e _ ai7r = (e r • e 01 ) (43) 

MP£ e _ ai7r = (e 7 * -QKcax-P) (44) 

At the ai vertex ai — > pir there are four relevant 4- vectors: e ai , e p , P ai , and Q ai = Pp — P^. 
These are formed into two independent amplitudes: 

= (e ai ■ e p ) (45) 

MEil^ = (e ai -Q ai )(e p -P ai ) (46) 

At the p vertex, there is only one matrix element 
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ME p ^ nn = (e p ■ Q p ) 



(47) 



as seen before. 

Combining all of these yields four independent amplitudes for the ayn process. These 



are 



ME Ln = 2 ( e 7* • e «i)( e ai ' e p)( e P ' Qp) 



(48) 



ME 1 ^ = ]T (ey- " • Q«i)(e p • P 0l )(e p • Q p ) 



(49) 



ME ai7T = ^ (e r ■ Q)(e ai • P)(e ai • e p ){e p ■ Q p ) 



(50) 



MP a ^ = £ • Q)(e ai ■ P)(e ai • Q ai )(e p • P ai )(e, • Qp) 



(51) 



Now, performing the sums over the polarizations yields: 



(52) 



ME 



ii 

<2l7T 



— (e 7 * • Q ai ) + 



(e r • P ai )(P i • Q 6 



p2 



~{Pa, ■ Q P ) + 



(P ai .P p )(P p .Q p )i 



Pp 2 



(53) 



MP a ^ = ]T(e T *.Q)[-(P.e p ) + 



(P-P ai )(P ai -e p ) 1 



P 2 

ai 



v^P • Qp) 



(54) 



MEll = (e r ■ Q) \-{P -Q ai ) + ( p - p ^-Q a .) 



<P ai -Q P ) + { - P ^ 



(P ai -P P )(P P -Q P ) 

>2 

P 



(55) 



Performing the sums over the polarizations of the p deferred in two of the four yields: 



MP, 



i f) \ , (e 1 *-Pp)(Pp-Q P ) 
— l e 7* ' VpJ H p5 



+ 



■P ai 



P 2 



-(Pa 1 -Qp)+ (PaVPp) pi PpQp) 

p 



(56) 
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ME. 



ii 

a\TT 



£7* ' Qai) + 



(^7* ' Pai ) (Pai ' Q 



ai , 



p2 

ai 



~{Pa, ■ Qo) + 



{P ai -P p ){P p -Q P h 



me!" = (t. 



aiir 



Q) 



-(P-Q P ) + {P - Pp) ^ Qp) 



+ 



p-p ai 

p2 



-{Pa, ■ Q P ) + {J ^0 



(Pa, -P P )(Pp-Qp) 
)2 

P 



ME™=( €r .Q) -(P.Q ai ) + 



{P-P ai ){P ai -Q ai ) 



PI 



■{Pa, ■ Qp) + 



{P ai -P p ){P p -Q p ) 



These matrix elements must be multiplied by Breit-Wigner formulas for a\ and / 
Bose symmetrized. 

* ai7r (Q) = R ie ^ (M< (nffi) + ME J ai MA)) 
+ R 2 eHMEi[{^A) + MEHMA)) 
+ R,eHMEl»{^l) + MEi[i{^)) 
+ R 4 e<t»{MEiY{nW) + ME^A)) 

E. The e+e" uj°tt° Mode 
This formalism can also provide a formula for the process 



e+e- - o/tt 



^0 ^ ^-^0 



The e + e — > c<j7r vertex has the matrix element: 
where P = P w + P^o and Q = P w — P^o 

The u; — > 7r + 7r~7r° vertex is governed by the matrix element. 



MP - e ^ P u P p P a 

iv± J-imrir ^pvpcr^u) 1 n+ tt~ u> 
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Since P w = P w + + P w - + P n o, this formula is equivalent to 

ME wnn = e, upa e»P» + P^P° (63) 

The complete matrix element is 

ME^ = Ea„ e^e^PPQ- (64) 
e a!3l 5€*P n+ P2-Pu (m 2 -m^-i^r^) ( 65 ) 

This formula reduces to a simple formula in terms of dot products of 4-vectors (polar- 
ization 4- vectors and 4- momenta). 

First, use the completeness relation for the uo to get 

ME^ = e^P^Q u e£e a ^ 5 P%PfPf 

(-g*a + (ml-ml-im^) ( 66 ) 

Note that: 



(VfhSpupn+pn pup 



P 



M 2 

since P w is repeated. Thus the matrix element becomes the simpler form: 



°def ~ < 



(67) 



ME^ = e^P^Q^t^P^P^Pf 

At this point, use the relation 

e^e^s = e^e pf3lS = (69) 

where 



if index repeated in abc or def 

1 if def is even permutation of abc (70) 
— 1 if def is odd permutation of abc 
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This is the Kronecker 3-delta. It is a generalization of the Kronecker delta Sj. 

1 if a = b 
otherwise 



(71) 



srab 
°cd 



The Kronecker 2-delta is defined by 

1 if ab is even permutation of cd 
— 1 if ab is odd permutation of cd 
if index repeated in ab or cd 
The Kronecker 2-delta can be expressed in terms of the Kronecker delta: 

C = (S a J b d - 5 a d 5 b c ) 



(72) 



(73) 



The Kronecker 3-delta can be expanded in terms of the simple Kronecker delta as well: 



m = mm - «/ - w c - w + w + w) 

This leads to an expansion in terms of dot products: 

A^Cptt^"^ 7 = ( A ■ D) (B • E) (C ■ F) - (B ■ D) (E ■ A) (C ■ F) 

- {A ■ D)(C ■ E){B -F)-(C- D)(B ■ E)(A ■ F) 
+ (B ■ D)(C ■ E)(A . F) + (C • D)(A ■ E){B ■ F) 



(74) 



(75) 



This can be used to produce a form for the un matrix element in terms of dot products 
of the 4-momenta and the virtual photon polarization 4-vector. 

ME = e r -[(P-P n -)(Q-P LU )P 7T+ 

- (P ■ P„){Q ■ P«-)P^ 

- (P ■ P« + ){Q ■ P„)P n - 
+ (P UJ -P)(Q-P 7r+ )P W - 
+ {P-P„ + ){Q-P„-)P UJ 

- (P ■ P n -)(Q ■ P W+ )P W ] 



(ml-m^-imuTu) 



(76) 
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V. CONCLUSION 



Some general formulas for simplifying Lorentz invariant amplitude expressions were de- 
veloped. These were applied to calculate the Lorentz invariant amplitudes for processes con- 
tributing to e + e~ — > 7r + 7r~7r°7r°. These formula are also applicable to e + e~ — > 7r + 7r~7r + 7r~ , 
e+e" -> K + K-tt°tt°, e+e~ -> K + K~K°K°, and e + e~ -> K + K~K + K~ . These formulas 
may be used to determine the electromagnetic form factors of the p meson and properties 
of the ai resonance. 
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TABLES 

TABLE I. Isobar Model Modes 
Mode 
e + e _ — > lott 
e + e~ — ► 
e + e~ — > ai7r 
e + e~ — > 7r + 7T _ 7r 7r° 
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